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Abstract. - We study dynamics of a classical particle in a one-dimensional potential, which 
is composed of two periodic components, that are time-independent, have equal amplitudes 
and periodicities. One of them is externally driven by a random force and thus performs a 
diffusive-type motion with respect to the other. We demonstrate that here, under certain 
conditions, the particle may move unidirectionally with a constant velocity, despite the fact 
that the random force averages out to zero. We show that the physical mechanism underlying 
such a phenomenon resembles the work of an escapement-type device in watches; upon reaching 
certain level, random fluctuations exercise a locking function creating the points of irreversibility 
in particle's trajectories such that the particle gets uncompensated displacements. Repeated 
(randomly) in each cycle, this process ultimately results in a random ballistic-type motion. 
In the overdamped limit, we work out simple analytical estimates for the particle's terminal 
velocity. Our analytical results are in a very good agreement with the Monte Carlo data. 



Introduction. - A long standing challenge has been a problem of how to rectify unbiased 
random or periodic fluctuations into directional motion. For microscopic fluctuations, induced, 
e.g., by Brownian noise, this question has been debated since the inception of the random walk 
theory, and several important concepts have been worked out PG] (see also Ref. Indeed, 
it has been established already by Curie 0] more than a hundred years ago that although 
a violation of the x — > —x symmetry is not sufficient to cause a net directional transport 
of a particle subject to a spatially asymmetric but on large scale homogeneous potential, the 
additional breaking of time reversal t — > — t symmetry, (e.g., due to dissipative processes), may 
lead to a macroscopic net velocity, such that here directed motion can result in the absence 
of any external net force. These systems, known currently as thermal ratchets [3], have been 
subject of intensive research, both theoretical (sec, e.g. Refs. |5I6I7I8I9I10I11I12I13I14I15I16| ) 
and experimental (see, e.g. Refs. ^I^IEIED!)- Advancements in this field have been 
summarized in Refs. |2Tll22[l23| . 

On the macroscopic scale, the same problem of obtaining a useful work from random or 
periodic perturbations has been known for many centuries and this task has been indeed 
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accomplished in many instances; to name but a few, we mention water- or windmills, or 
watches. In the latter case, the watch makers had to create a mechanism capable to convert 
the raw power of the driving force into regular and uniform impulses, which was realized by 
inventing various so-called escapement devices (see, e.g., Ref. |24| ) . Such a devise is most 
often a sort of a shaft or an arm carrying two tongues - pallets, which alternately engage 
with the teeth of a crown-wheel. The pallets follow the oscillating motion of the controller - 
the balance-wheel or a pendulum, and in each cycle of the controller the crown-wheel turns 
freely only when both pallets are out of contact with it. Upon contacts, the pallets provide 
impulses to the crown-wheel, (which is necessary to keep the controller from drifting to a 
halt), and moreover, perform a locking function stopping the train of wheels until the swing 
of the controller brings round the next period of release. 

In this Letter we exploite the concept of such an escapement device to build a molecular 
motor capable of doing a " useful work" under the action of random forces which average out to 
zero. The system we consider consists of a classical particle in a one-dimensional (ID) two- wave 
potential, which is is composed of two periodic in space, identic time- independent components. 
One of them is externally driven by a random force and thus performs a random translational 
motion with respect to the other. This model has been previously proposed in Ref. |14j and it 
has been discovered that under certain conditions and for constant or periodic in time external 
driving force, the particle may perform totally directed motion with a constant velocity. Here 
we demonstrate that in such a system the particle may move unidirectionally with a constant 
velocity even in situations when the direction of the driving force (the controller) fluctuates 
randomly in time. We show that the physical mechanism underlying such a phenomenon is 
indeed that of an escapement-type device - upon reaching certain levels, random fluctuations 
exercise a locking function creating the points of irreversibility in particle's trajectories such 
that, (despite the fact that the fluctuations average to zero), the particle gets uncompensated 
displacements. Repeated (randomly) in each cycle, this process ultimately results in a ballistic- 
type motion. Focussing on the overdamped limit, we use this physical picture to map the 
original system onto a Brownian motion (BM) process on a hierarchy of disconnected intervals. 
This allows us to work out simple analytical estimates for the particle's terminal velocity, which 
are in a very good agreement with Monte Carlo data. 

The Model. - Following Ref. |14j . we consider a simple piece- wise continuous potential 



where the parameter £ £ (0,1) determines the asymmetry of the potential, with £ = 1/2 
corresponding to the symmetric case. We will constrain our analysis here to the case £ < 1/2. 
Note that it has a periodicity b = 1, so that H(x + 1) = Tl(x) \/x, an amplitude n = 
maxn(x) = — minll(:r), and one minimum is located at x = 0, i.e. 11(0) = — n . 

Now, the total potential V(x, 7) is the sum V(x, 7) = H(x) + H(x — 7), where 7 defines the 
external translation. Because of the periodicity of the potentials n(x), the potential V(x, 7) 
is periodic in both arguments, so that V(x, 7 + 1) — V(x, 7) V7. In this potential landscape, 
the equation of motion of a particle of mass m in a ID medium with damping n reads: 




(1) 



mx t + nx t + 



dV{x ul ) 
dx t 



= 0, 



(2) 
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Xt being the particle's trajectory. We note that energy is steadily pumped into the system 
through 7, and is dissipated when r\ > 0, which prevents particle's detachement. 

Finally, we specify the properties of 7. We suppose that 7 obeys the Langevin equation, 
7 = ft, where ft is a random Gaussian, delta-correlated force with moments 

Tt = fo, T~h=f 2 +2D8(t-t'), (3) 

while the overbar denotes averaging over thermal histories. We will focus here mostly on the 
case /o = 0, such that 7 is just a trajectory of a symmetric ID BM. The case when /o ^ or 
when 7 represents anomalous diffusion will be briefly discussed at the end of the paper. 

Particles Dynamics in the Overdamped Limit. - We will restrict ourselves to the limit 
of an overdamped motion, rj/[(2Tr) -y/milo] 3> 1. As shown in Ref. |14| . in this limit dynamics 
of the particle is governed entirely by the evolution of different minima of the total potential 
V(x, 7). Consequently, in order to obtain the trajectory xt and the average velocity V = it, 
it suffices to study the time evolution of positions of these minima. This process has been 
amply discussed in Ref. [14] and here we will merely outline the main conclusions. 

According to Ref. 14 , the particle's dynamics proceeds as follows: the total potential 
V(x, 7) possesses a set of minima and position of each minima changes in time as the transla- 
tion 7 evolves. The particle, located at t — at the first minimum simply follows the motion 
of this minimum up to a certain moment of time, or more precisely, a time moment when 
this minimum reaches a certain point x = x. These appoints are the points of instability X 
in the (x, 7) plane, where the corresponding local minimum of the potential V(x, 7) ceases 
to exist, and which emerge due to the asymmetry of IL(x). Further on, when the minimum 
disappears, a particle located at such a point performs an irreversible motion jumping to one 
of two neighboring minima, which still exist. Since the potential in the vicinity of each x may 
have both a downward and an upward bending, the jumps may be performed in both the 
leftward and the rightward directions. Depending whether £ < 1/2 or £ > 1/2, the leftward 
or rightward jumps may result in uncompensated or compensated displacements. 

Evolution of the effective translation. - We now put the particle's dynamic rules deduced 
in Ref. |14j into a mathematical framework, amenable to analytical analysis, and introduce an 
auxiliary stochastic process 74. We will call 74 as an "effective translation", since it follows, 
apart of some special points (to be specified below), the evolution of the potential minima 
visited by the particle and keeps track of both the translation 7 and the effect of the irre- 
versibility points which induce particle's leftward and rightward jumps to neighboring minima 
when the local minimum seizes to exist. We note also that the particle's actual trajectory 
xt = Jt almost everywhere, except for some special points (to be specified below). 

The effective translation 74 obeys the following Smoluchowski-Feynman ratchet form: 

^=t{ lt ) + f t , (4) 

where f t are random forces defined in cq.Q, while C, (74) is a local operator, which may be 
thought of as a space-derivative of some "effective" potential. 

This operator equals zero everywhere except for the set of special points jt = ~K + £ and 

7t = —K + 3/2, K = 1, 2, At these special points, the action of the operator C (74) is as 

follows: When j t reaches jt = —K + 3/2, (which we call the "reflection point"), the operator 
C (7t) changes 7* — > 7t — 1, (i.e. shifts its value from 7 t = —K + 3/2 to j t = —K + 1/2). 
On the other hand, when j t hits jt = —K + £, the operator C (74) shifts the position of the 
reflection point from — K + 3/2 to —K + 1/2, thus performing a "locking" function - 74 can 
not now return below the point — K + 3/2. 
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Fig. 1 - Typical realization of the effective translation 7*. The dashed lines denote the "reflection 
points" x — —K + 3/2 of the trajectory jt, while the dash-dotted lines - the points — K + £. In 
simulations, we took £ = 1/4, the lattice spacing £7 = 2 x 1CF 2 and the characteristic jump time 
t = 1; hence, D = (Sjf /2r = 2 x 1(T 4 . 

Fig. 2 - Time evolution of 74 viewed as a BM on a hierarchy of disconnected intervals. 



In Fig.l we depict a typical realization of the process 74 for /o = 0. Note that because of 
C ("ft), It experiences an effective drift in the negative direction (for £ < 1/2); for C (jt) = 
one should, of course, regain a standard result 77 = 7 = and = 7 2 = 2Dt. 

Note now that j t can be viewed from a different perspective; namely, 74 can be regarded 
as a BM on a hierarchical lattice composed of a semi-infinite set of disconnected intervals 
of length 3/2 — £ (see Fig. 2). The process starts at t = at the origin of the first interval 
(K = 1) and evolves freely until it cither hits the right-hand-side (RHS) boundary j t = 1/2 
(the reflection point), in which case it gets transferred instantaneously at position 74 = —1/2 
and continues its motion on the interval K = 1, or reaches the left-hand-side (LHS) boundary 
7t = — 1 + £ - the "exit point" and gets irreversibly transferred to the second (K = 2) interval. 
In the second, and etc, interval 74 evolves according to the same rules. 

The particle's trajectory x t follows the process 74, i.e. x t = Jt, except for two details: 
a) when 7 t e [-K + 1,-K + 3/2], K = 1,2,..., (thick lines in Fig. 2), x t sticks to the 
left boundary of these intervals, such that the trajectories xt are saltatory with the random 
pausing times corresponding to the times spent by j t in the intervals [— K + 1, — K + 3/2], 
K = 1, 2, . . ., (see Fig.l). b) the event when j t hits the LHS boundary of the K-th interval 
and gets instantaneously transferred to the next interval (i.e. 7* does not have a discontinuity) 
corresponds to the event when the particle makes a jump from x t = —K + t; to x t = —K, i.e. 
at this point x t changes discontinuously. On the other hand, it is clear that such jumps do 
not contribute to the average particle velocity x t and hence, x t =V = / y t - 

Particle's terminal velocity. - We turn next to the computation V for / = 0. Evidently, 
even in this simplest case the map depicted in Fig. 2 is too complex to be solved exactly. 
However, some simple analytical arguments can be proposed to obtain very accurate estimates 
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of V. To do this, we first note that it is only when passing from the K-th interval to the 
(K + l)-th one, (which is displaced on a unit distance on the 7t-axis with respect to the 
previous one), the process 7< gains an uncompensated (negative) displacement equal to the 
overlap distance of two consequitive intervals, i.e., 1/2 — £. Consequently, we may estimate V 
as V = —(1/2 — £)/T, where T is the mean time which 74 "spends" within a given interval. 
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Fig. 3 - Time evolution of the effective external translation 7 t on the K-th interval. 

Fig. 4 - The absolute value of V versus £. Solid line gives the analytical prediction in eq.ljjjj, while 
the triangles depict Monte Carlo data. Note that we use, for notational convenience, the log-log scale. 



To estimate T, let us consider how 74 evolves on a given interval; T is, as a matter of fact, 
the mean time needed for j t to reach diffusively for the first time the LHS boundary of the 
interval - the "exit" point —K + £, starting from the "entry" point —K + 1 + £ (see Fig. 3). 
We divide next the K-th interval into three sub-intervals: l\ = I3 = 1/2 — £ and I2 = 1/2 + £, 
Fig. 3. Since l\ = I3 < I2, (£ > 0), for typical realizations of ft, the following scenario should 
hold: 74 evolves randomly around the "entry" point and first hits the RHS boundary of the 
K-th interval passing thus through the sub- interval l\; the mean time required for the first 
passage of h is denoted as T x . For BM with diffusion coefficient D, T x = (1/2 - £) 2 [ /2D 
Further on, after hitting the RHS boundary, j t gets reflected to the point jt = —(2K — l)/2 
and evolves randomly around this point until it hits the LHS boundary of the interval - the 
"exit" point. For symmetric BM, the time X3 needed for the first passage through I3 is equal 
to T\. Consequently, T = Ti + T 3 = (1/2 - £) 2 /D and the velocity reads 

V = -2D {I - 2£) _1 (5) 

In Fig. 4 we compare our prediction in eq.JHJ and the results of Monte Carlo simulations, which 
shows that theoretical arguments presented here capture the essential physics underlying the 
dynamics of the random map depicted in Fig. 2 and, hence, of the process in eqs.© and |J3). 
Now, several comments on the result in eq.Q are in order. 

(i) First, we note that V, eq.lJSJ, diverges when £ — > 1/2 (note, however, that V = for 
£ = 1/2), which is a seemingly counter-intuitive behavior. Such a behavior stems, as a matter 
of fact, from our definition of j t ', namely, by writing our eq.Q, we tacitly assume that both 
the effective step J7 of the process jt and the characteristic jump time r are infinitesimal 
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variables, (while the ratio D — (5j) 2 /2t is supposed to be fixed and finite). On the other 
hand, for any real system $7 and r might be very small but nonetheless are both finite. For 
finite Sj and r, j t is a symmetric hopping process on a lattice with spacing £7 and with the 

characteristic time r. Then, we have Ti = T 3 = r(l + LjL, where L = (1/2 - £)/Sj is the 

number of elementary steps £7 in the interval (1/2 — £). This yields, in place of eq.JSJ, 

V = -2D (28j + 1 - 2^,)~ 1 , (6) 

where now V tends to a finite value when £ — > 1/2 (but still V = when £ = 1/2). 

(ii) Second, notice that V in eq.JSJ) is a monotoneously increasing function of £; that is, 
the absolute value of V is maximal when £ — > 1/2 and minimal, |V| = 2D, for the strongest 
asymmetry, £ = 0. We note that such a behavior again stems from the definition of ft as a 
Gaussian, delta- correlated noise. On the other hand, since ft is influenced by some external 
processes, it might be characterized, e.g., by correlations or be a non-local in space or in time 
(discontinuous) stochastic process. In both cases, the behavior of V in the presence of such 
an external force would be different of that predicted by eqs.JSJ) and ©. 

Consider now what will happen if f t represents the so-called delta-correlated Levy process 
(LP) with parameter /j, (see an exposition in excellent Ref. [23] for more details). We note 
here parenthetically that the case /i = 2 corresponds to the usual Gaussian case, which yields 
conventional BM, while the case [i = \ describes the so-called Cauchy process. 

Now, what basically changes when we assume that f t is the LP, is the form of the first 
passage times T\ and T 2 . Here, the time required for the first passage of an interval of length 
1/2 - £ reads T x = T 3 ~ (1/2 - 25 and hence, V ~ -(1 - 2£) 1 -' J . 

Consequently, we infer that V will be a monotoneously increasing function of the parameter 
£ only for the LP with fj, > 1, i.e. for the persistent processes (which have a finite probability 
for return to the origin). In the bordeline case of /1 = 1 velocity will be independent of the 
asymmetry parameter. On contrary, for transient (having zero probability of return to the 
origin) processes with fj, < 1, which are not space-filling (fractal) and occupy the space in 
clustered or localized patches, one would find that V is a decreasing function of £. 

(iii) Finally, we consider the case when the bias is oriented in the positive direction, i.e. 
fo > 0. In this case, 74 is an asymmetric hopping process on a ID discrete lattice of spacing 
$7 and with transition probabilities p ( q ) of jumps in the positive (negative) directions which 
obey p/q = exp(/3/ ^7)j with p + q = I, where (3 denotes the reciprocal temperature. 

Here, the symmetry T\ — T3, which exists for fo — 0, is broken: when passing through l\ 
the process 74 follows the field, while the passage through I3 takes place against it. Supposing 
that ^7 and r are both finite, we have that here: 



Ti,3 = 



p-q 




(7) 



where the upper (lower) sign corresponds to index "1" ("3"), while (pi — q and ^3 = p. Note 
that T\ grows linearly with the sub- interval length 1/2 — £, while T3 shows much stronger, 
exponential interval-length dependence, and hence, controls the overall time spent within a 
given interval. Consequently, the particle's velocity in this case attains the form 



(1/2-0 sinh^(/3/ 57/2 

V = 



cosh (W7/2) [cosh (/3/ (l - 2£ + 6 7 )/2j - cosh 



(8) 
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which reduces in the diffusion limit to the following result 



V = - 



(1 - 2£) J D/? 2 /o 2 



(9) 



8sinh 2 (/3/ (l - 20/4 



) 



The salient feature of these results is that here the particle's drift proceeds against the applied 
field, due to fluctuations in the process jt - some (exponentially small) number of trajectories 
which travel against the field. 

In conclusion, we have studied dynamics of a classical particle in a ID potential, composed 
of two periodic components, one of which is driven by an external random force. We have 
shown that in such a system the particle may move unidirectionally with a constant velocity 
even when the random driving force averages out to zero. We have demonstrated that the 
physical mechanism underlying such a behavior resembles the work of the so-called escape- 
ment device, used by watch makers to convert the raw power of the driving force into uniform 
impulses; here, indeed, upon reaching certain levels, random forces lock the particle's motion 
creating the points of irreversibility, such that the particle gets uncompensated displacements. 
Repeated (randomly) in each cycle, this process ultimately results in a ballistic-type motion. 
Concentrating on the overdamped limit, we have used this picture to map the original system 
onto a BM process on a hierarchy of disconnected intervals, which allowed us to present ana- 
lytical estimates for the particle's velocity. Our analytical results are in a very good agreement 
with Monte Carlo data. Extensions for systems with other than Gaussian fluctuations and 
systems with global bias have also been presented. 
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